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We study the thermodynamics of the one-dimensional extended Hubbard model at half-filling 
using a density-matrix renormalization group method applied to transfer matrices. We show that 
the various phase transitions in this system can be detected by measuring usual thermodynamic 
quantities like the isothermal compressibility and the uniform magnetic susceptibility. For the 
isothermal compressibility we show that universal crossing points exist which allow to accurately 
determine the line where the charge gap vanishes. By studying in addition several correlation 
functions, we confirm the existence of a phase with long-range dimer order (bond order) which has 
been a matter of debate for several years. According to our calculations this phase is located in 
a narrow region between the spin-density and charge-density wave phases up to a tricritical point 
which we estimate to be at Ut = 6.7 ± 0.2, Vt = 3.5 ± 0.1. Our results for the phase diagram 
are in good agreement with the most recent zero-temperature density-matrix renormalization group 
study, however, they disagree in some important aspects from the most recent Quantum-Monte- 
Carlo study. 
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PACS numbers: 71.10.Fd, 71.10.Pm, 71.10.Hf, 05.70.-a 



I. INTRODUCTION 

Understanding the effects of competing interactions 
and the associated quantum phase transitions in mod- 
els of strongly correlated electron systems is still one of 
the main issues of modern condensed matter physics. In 
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one dimension correlation effects are particularly strong 
and a number of analytical and numerical tools have been 
developed for this case. This has led to intense work on 
such systems in the last decades. One of the seminal 
models in this context is the extended Hubbard model 
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Here Cj ^ {cj^a) is a creation (annihilation) operator for an 

electron with spin cr | at site j, Uj^cr = '^ja'^3,<y^ ^^^"^ 
rij = Tij^i + rij^i. t is the amplitude of a nearest neighbor 
hopping, C/ > an on-site and > a nearest-neighbor 
Coulomb repulsion. A possible additional magnetic field 
is denoted by h and the chemical potential by fi. Here 
we will concentrate on the half-filled case n — 0. In 
the following, we will only keep the hopping amplitude t 
where it is necessary for the sake of clarity and set t = 1 
otherwise. 

In the strong coupling limit, U,V ^ t, it is easy to 
see by simple energetical considerations that two differ- 
ent ground states exist: For U < 2V the system is an 
insulator with long-range charge density wave (CDW) 
order whereas for U > 2V a state with quasi-long-range 
spin density wave (SDW) order forms. The transition be- 
tween these two phases in the strong coupling limit is first 
order. In the weak coupling limit, U,V <^ t, the model 
can be studied using bosonization and g-ology.^'^ In this 
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framework one finds again a phase transition between 
the SDW and CDW phase at U = 2V. In the spin sector 
this transition is driven by an operator which turns from 
marginally irrelevant in the SDW phase to marginally rel- 
evant in the CDW phase. The spin gap therefore opens 
up exponentially slowly and the transition in the spin 
sector is of Kosterlitz-Thouless (KT) type. In the charge 
sector, on the other hand, there is a relevant operator in 
both phases leading to a charge gap. The amplitude of 
this operator vanishes only at the transition line U — 2V 
so that the charge gap disappears.* The transition in the 
charge sector is therefore second order. Already from 
the strong and weak coupling approaches it is clear that 
a point {Ut, Vt) in the intermediate coupling regime must 
exist where the order of the phase transition changes. 

In the last few years the extended Hubbard model 
has attracted renewed attention because it has been sug- 
gested that the phase diagram obtained by the weak cou- 
pling g-ology approach and strong-coupling perturbation 
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theory might not be complete. Nakamura pointed out 
first that there is no symmetry requiring the hnes in 
C/,y-parameter space, where tlic marginal operator in 
the spin sector changes sign and where the relevant oper- 
ator in the charge sector vanishes, to coincide. ^'^ The cou- 
pling constants for these two operators do coincide in the 
standard g-ology approach where they are calculated to 
first order in the interaction parameters. However, they 
might differ once higher order corrections are taken into 
account. This opens up the possibility for an intermedi- 
ate phase. By extracting the scaling dimensions related 
to the critical exponents of certain correlation functions 
from finite size energy spectra, Nakamura indeed found 
a phase with long-range dimer order in a small region 
between the SDW and CDW phases. This phase is often 
called a bond-order wave (BOW) state. The existence 
of such a phase around U = 2V in the weak coupling 
regime was supported by Quantmu-Monte-Carlo (QMC) 
calculations'"'* as well as by a g-ology approach where the 
coupling constants have been calculated beyond leading 
order. ^ However, in a first DMRG calculation^" such a 
phase was only found above the tricritical point (Ut, Vt) 
and only directly at the first order transition line. A 
later DMRG calculation,^^ on the other hand, qualita- 
tively confirms again the phase diagram as proposed by 
Nakamura. Further evidence for the existence of a BOW 
phase in the weak coupling regime was also provided by 
a functional renormalization group analysis. 

Although the most recent DMRG"'^"'^ and the most re- 
cent QMC study^ agree that a BOW phase of finite ex- 
tent does exist, they disagree about the shape of this 
phase. Whereas in the phase diagram of Ref. 11 the BOW 
phase ends at the tricritical point, it extends beyond this 
point to larger values of U, V in the phase diagram of 
Ref. 8. The question whether the tricritical point also 
marks the end of the BOW phase or is located on the 
BOW-CDW boundary therefore remains an open issue. 

In this work we will investigate the half-filled one- 
dimensional extended Hubbard model using a density- 
matrix renormalization algorithm applied to transfer ma- 
trices (TMRG) . This numerical method allows it to calcu- 
late thermodynamic properties of the model in the ther- 
modynamic limit. We will provide further evidence for 
the correctness of the phase diagram as first proposed by 
Nakamura and give an estimate for the tricritical point 
(f/(,Vt). In particular, we will argue based on our nu- 
merical results that the BOW phase ends at the tricrit- 
ical point and does not extend to larger values of U, V 
in contrast to the findings in Ref. 8. In the process, we 
will develop and discuss criteria to identify the different 
phases and transition lines by considering usual thermo- 
dynamic quantities like the uniform magnetic susceptibil- 
ity, the isothermal compressibility (charge susceptibility), 
and the specific heat. 

Our paper is organized as follows: In Sec. II we briefly 
introduce the TMRG algorithm and compare results for 

the Hubbard model {V = 0) with exact results obtained 
by the Bethe ansatz. In Sec. Ill we then present results 



for a variety of thermodynamic quantities which allow us 
to determine the phase diagram of the extended Hubbard 
model at half-filling. The last section is devoted to our 
conclusions. 



II. THE TMRG ALGORITHM AND THE 
HUBBARD MODEL 

The density-matrix renormalization group applied to 
transfer matrices (TMRG) has been explained in detail in 
[13, 14, 15]. Here we only want to briefly discuss the most 
important aspects. The TMRG algorithm is based on a 
mapping of a one-dimensional quantum system to a two- 
dimensional classical one by means of a Trotter-Suzuki 
decomposition. In the classical model one direction is 
spatial whereas the other corresponds to the inverse tem- 
perature. For the classical system a so called quantum 
transfer matrix (QTM) is defined which evolves along the 
spatial direction. At any non-zero temperature the QTM 
has the crucial property that its largest eigenvalue Aq 
is separated from the other eigenvalues by a finite gap. 
The partition function of the system in the thermody- 
namic limit is therefore determined by Aq only, allowing 
it to perform this limit exactly. The Trotter-Suzuki de- 
composition is discrete so that the transfer matrix has a 
finite number of sites or local Boltzmami weights M. The 
tempcaature is given by T {eM)~^ where e is the dis- 
cretization parameter used in the Trotter-Suzuki decom- 
position. The algorithm starts at some high-temperature 
value where M is so small that the QTM can be diagonal- 
ized exactly. Using a standard infinite-size DMRG algo- 
rithm, sites are then added to the QTM leading to a suc- 
cessive lowering of the temperature. A source for a sys- 
tematic error in these calculations is the finite discretiza- 
tion parameter e. However, this only leads to errors of or- 
der in all thermodynamic quantities considered in the 
following. We will choose e = 0.025 or 0.05 so that this 
systematic error will only be of order 10^"^ — 10^^. An- 
other error is introduced by the truncation of the Hilbert 
space. This error will grow with decreasing temperature 
and will finally make the calculations unreliable. Down 
to which temperature the DMRG algorithm works will 
depend on the maximum number of states N kept in the 
truncated Hilbert space basis. The truncation error is dif- 
ficult to estimate. We therefore start by comparing our 
TMRG results for the Hubbard model (V = 0) with exact 
results obtained by Bethe ansatz. Within the TMRG 
algorithm nothing changes fundamentally when we intro- 
duce the nearest-neighbor Coulomb repulsion V so that 
we expect a similar accuracy in this case. 

As an example, we consider the case {7 = 8. Results 
with a similar accuracy are also obtained for other U. 
Using the TMRG method, the free energy per site is given 

by 

/ = -TlnAo. (2) 

The specific heat is then obtained by C = -Td'^f/dT'^ 



FIG. 1: (Color online) TMRG data for the specific heat C 
of the Ilubbard model a,t U — 8 with = 200 states kept 
(red solid line) compared to Bethe ansatz data (circles) as a 
function of temperature T. The lower graph shows the error 
AC of the TMRG calculation. 
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FIG. 2: (Color online) Same as Fig. 1 for the magnetic sus- 
ceptibility Xs- 



and is shown in Fig. 1. It is also easy to calculate the 
expectation values of local operators with the TMRG al- 
gorithm. To obtain the magnetic susceptibility, Xs: the 
expectation value m = (S*^) — (n| — n[)/2 is calculated 
in the presence of a small magnetic field 6h ~ 10~^. The 
susceptibility is then given by Xs = m/Sh and shown in 
comparison to the exact result in Fig. 2. Similarly, the 
isothermal compressibility (charge susceptibility), Xc, is 
obtained by applying a small chemical potential and 
is shown in Fig. 3. For the spin and charge susceptibil- 
ities Xs, Xc the error does not exceed 5 x 10"** down to 
temperatures T « 0.05. For the specific heat C, the er- 
rors are an order of magnitude larger because a second 
numerical derivative has to be calculated. 



III. THE PHASE DIAGRAM OF THE 
EXTENDED HUBBARD MODEL AT 
HALF-FILLING 

To investigate the phase diagram we will consider a 
number of different thermodynamic quantities like mag- 



FIG. 3: (Color online) Same as Fig. 1 for the charge suscep- 
tibility Xc- 



netic susceptibilities, compressibilities, specific heats, 
and expectation values of local operators. Furthermore, 
we will study the behavior of correlation lengths which 
can be obtained within the TMRG method by calculat- 
ing next-leading eigenvalues of the QTM.^^'^^ Depending 
on the required accuracy and the temperature regime we 
want to access the basis for the truncated Hilbert space 
wiU consist of TV = 200 - 400 states. 

We start with the strong coupling limit where the two 
existing phases and the first order phase transition be- 
tween these phases are well understood. We then derive 
an estimate for the tricritical point where the first or- 
der line ends. Next we discuss the considerably more 
complicated weak coupling regime and present the phase 
diagram as obtained by TMRG. Finally, we will address 
the controversial question whether or not the BOW phase 
ends at the tricritical point. Throughout, we will discuss 
in how far one can identify the different phases and phase 
transitions by studying only easily measurable thermody- 
namic quantities like the specific heat, magnetic suscep- 
tibility and compressibility. 



A. Strong coupling 

In the strong coupling limit, U,V 3> t, the ground 
state energy can be systematically expanded in terms of 
the hopping parameter t. In lowest order, the hopping 
can be completely neglected. Then, depending on the ra- 
tio U/V, two different ground states are possible. These 
states are depicted in Fig. 4. The energy of the CDW 
state is then given by E^jj^y ~ LU/A — LV with L being 
the number of lattice sites. The SDW state has energy 
^SDW ~ ~LU /A. The two energies as a function of [/, V 
cross at [/ = 2V resulting in a first order phase transition. 
As usual, in second order in t, virtual hopping processes 
lead to an effective antiferromagnetic coupling of Heisen- 
berg type for the spins in the SDW state with coupling 
constant J = 2t^/{U — V)}''^ This state therefore has 
a charge gap but no spin gap and algebraically decay- 
ing spin correlation functions. The CDW state, on the 
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CDW: (li)0(n)0(!i)0(!i)0 

SDW: ©©©(D®©®0 

FIG. 4: The two ground states in the strong couphng hmit 
U,V ^ t. The state in the first line is a CDW state where 
every second site is doubly occupied, whereas the state in the 
second line is a state with every site singly occupied. Virtual 
hopping processes induce a quasi long-range SDW order for 
this state. 
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FIG. 5: Left column: (1) CDW state with one particle added, 
(2) CDW state with one double occupancy broken. Right 
column: (1) SDW state with one particle added, (2) SDW 
state with one double occupancy. 



other hand, has a charge and a spin gap. Excitations for 
the CDW and SDW state, ignoring hopping processes, 
are shown in Fig. 5. In lowest order perturbation the- 
ory, the energies of the excited states depicted in Fig. 5 
are given by E^jjy^, = E^jy^ - U/2 + 2V, EIj^^^ = 
^CDW — U + iV for the excited CDW states, and 
EIdw = E'sDW + U/2, El^^ = El^^ + U -V iov the 
excited SDW states. Excitation (1) in the CDW state 
is a charge excitation, whereas the breaking of a double 
occupancy - excitation (2) - is a spin excitation. If we sep- 
arate the two single spins in this excited state we obtain 
an excitation energy E^^y^ + 2[—U/2 + 2V), i.e., each 
single spin contributes —U/2 + 2V. In thermodynamic 
data the activated behavior will be characterized by the 
energy of a single excitation irrespective of whether these 
excitations appear in pairs or not. In the strong coupling 
limit, it follows that at the transition line charge and spin 
gap as obtained from thermodynamic data are expected 
to be equal As = Ac = U/2 and that both gaps increase 
linearly ~ 2V away from the transition line. In the SDW 
phase, excitation (1) is also a charge excitation and has 
a lower energy than excitation (2). The charge gap in 
the SDW phase is therefore given by Ac = U/2 and is 
independent of V. 

In Fig. 6, TMRG results for the spin susceptibility Xs 
and the spin gap As at U = 12 are shown. If a gap 
A exists, the dispersion of the elementary excitations is 
given by e{k) ~ A-|-A;^/(2m) with some effective mass to. 
It is then easy to see that the corresponding susceptibility 
will show activated behavior 

x(T)~22t|ZI) ,3) 
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FIG. 6: (Color online) Magnetic susceptibilities (solid lines) 
for 17 = 12 and y = 1.0, 2.0, • • ■ , 6.0, 6.1, 6.2, • • • ,7.0 as a 
function of temperature T. The dashed lines are fits according 
to Eq. (3). The inset shows the spin gap As extracted from 
those fits (circles) as a function of V. The solid lines in the 
inset denote the theoretical result in the strong coupling limit. 
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FIG. 7: (Color online) Charge susceptibilities (circles) for U — 
12 and V = 6.0, 6.1, . . . , 7.0 as a function of temperature T. 
The lines are a guide to the eye. The dashed lines are fits 
according to Eq. (3). The inset shows the charge gap Ac 
extracted from those fits (circles) as a function of V. The 
solid lines in the inset denote the theoretical result in the 
strong coupling limit. 



at temperatures T <C A. Using this function to fit the 
numerical data, we are able to extract the spin gap As- 
As shown in the inset of Fig. 6 the behavior of As as 
a function of y at C/ = 12 is already reasonably well 
described by the strong coupling limit, i.e., there is no 
spin gap up to « U/2, then A^ jumps to approximately 
U/2 and then increases linearly with slope 2. 

Similarly, we show TMRG results for the charge sus- 
ceptibility Xc and the charge gap Ac at U — 12 in Fig. 7. 
The results obtained for the charge gap Ac are also al- 
ready close to the strong coupling limit, although the 
gap is a bit smaller than U/2 in the SDW phase and it 
shows some V dependence when the transition point is 
approached. 
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FIG. 8: (Color online) Double occupancy as a function of 
temperature for U — 12 and different V. Inset: Extrapo- 
lated values of the double occupancy at zero temperature as 
a function of V. 



FIG. 9: (Color online) Specific heat as a function of temper- 
ature for U = 12 and different V. Inset: A sharp peak forms 
just above phase transition. 



Another quantity which allows us to detect the phase 
transition and to determine its order is the double occu- 
pancy 

d = {njAnj^i) . (4) 

In the strong coupling limit at zero temperature d = in 
the SDW state and d=l/2 in the CDW state. In Fig. 8 
we show d for U = 12 and various V. In the extrapolated 
data for zero temperature some corrections to the strong 
coupling limit are visible, d is already nonzero in the 
SDW phase and increases slightly with V. However, a 
jump in d at y « 6.18 is obvious. In the CDW phase d 
continues to increase with V and approaches 1/2 in the 
large V limit. 

The specific heat shown in Fig. 9 has two maxima for 
U = 12 and V — 0. The lower and higher tempera- 
ture maximum are due to spin and (gapped) charge ex- 
citations, respectively.^^ At low temperatures only the 
gapless spin excitations do therefore contribute and con- 
formal field theory predicts 



With increasing V the spin velocity Vs increases leading 
to a decreasing slope and to a shift of the lower tempera- 
ture maximum to higher temperatures. At the same time 
the charge gap decreases leading to a shift of the higher 
temperature maximum to lower temperatures. The be- 
havior changes drastically above the phase transition 
V > Vc ~ 6.18, because in the CDW phase the spin exci- 
tations are now also gapped and the specific heat shows 
activated behavior C ^ e^'^l^ with A = min(As,Ac). 
The emergence of a sharp peak for V can be un- 

derstood as follows: Because C{T)dT = — eo with Co 
being the ground state energy, the area under the curve 
will be nearly unchanged when going from a value just 
below the phase transition, say I^ = 6.1, to a value just 



above the transition, say V — 6.2. In addition, also the 
high temperature behavior will be almost unaffected by 
this small change in V. Hence the weight suppressed by 
the gap at low temperatures will show up in a sharp peak 
just above the gap. This is shown in the inset of Fig. 9 
and constitutes one possibility to detect the first order 
transition easily from thermodynamic data. 

B. The tricritical point 

From the discussion in the introduction it is clear that 
the first order transition line must end at some point 
{Ut,Vt) because the phase transitions at weaker couplings 
are expected to be continuous. We found that a good 
criterion to determine this endpoint with the TMRG 
method is to study the double occupancy d (4) . As shown 
for the case [/ = 12 in Fig. 8, d as a function of T shows 
a dramatically different behavior depending on weather 
we choose a V such that we are in the SDW phase or 
a V such that we are in the CDW phase, d at a fixed 
U extrapolated to zero temperature therefore shows a 
jump Ad as a function of V if the phase transition is first 
order. Reducing the on-site repulsion U we expect this 
jump to become smaller and smaller until it disappears at 
Ut- For U = 7.0 we can still detect a finite jump ~ 0.17 
at F ~ 3.65 (see Fig. 10(c)) whereas d as a function of 
V seems to be continuous for t/ = 6 (see Fig. 10(b)). 
To determine the point ([/< , Vf ) more accurately we have 
plotted the jump A^ as a function of U in Fig. 10(d). 
We can fit these data very well by a power law which 
leads us to the estimate Ut = 6.7 ± 0.2. Because the 
value for [/ — 7 is least reliable, we also did fits where 
this point was excluded. Similarly, we tried fits where 
the data points for large U were excluded. The results 
of the various fits lead to the error estimate above. For 
each possible value of Ut we can find Vt with high ac- 
curacy. For the values of Ut estimated above, we have 
Vt = 3.5 ±0.1. Here the uncertainty in Vt is not an error 
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FIG. 10; (Color online) Extrapolated values of the double oc- 
cupancy d at zero temperature as a function of V for different 
U. The inset of (c) shows a zoom of the region where d jumps, 
(d) Extrapolated TMRG data (circles) for the jump Ad in the 
double occupancy at the phase transition at zero temperature 
as a function of U. The line is a fit = 0.232(t/ - 6.7f''^^. 



estimate but rather means that Vt ~ 3.4 for Ut — 6.5 and 
Vt « 3.6 for Ut = 6.9. 



C. Weak coupling 

The phase diagram in the weak coupling limit is more 
complicated than in the strong coupling limit. Instead of 
a first order, we expect different continuous phase tran- 
sitions here. Theoretically, the weak-coupling limit can 
be investigated by bosonization with the coupling con- 
stants of the operators in the effective Hamiltonian being 
determined in first order in the interaction parameters. 
This method is also often termed g-ology.^ As usual in 
one-dimension, the charge and the spin sector completely 
separate in the low-energy effective bosonic theory due 
to the linearization of the excitation spectrum. In the 
charge sector at half-filling, Umklapp scattering leads to 
a relevant interaction term in the bosonic Hamiltonian 
which creates a charge gap. In the spin sector, on the 
other hand, the leading interaction term corresponding 
to backward scattering is only marginal. The amplitudes 
of both terms in the weak-coupling limit are proportional 
to U—2V.'*' The system therefore has always a charge gap 
except at U — 2V where the amplitude of the Umklapp 
scattering term vanishes. The charge gap at fixed U near 
the phase transition behaves as 



A. 



W-VA 



(6) 



with a > being an interaction dependent critical expo- 
nent and Vc ~ U/2 at weak coupling. This means that 



the transition in the charge sector is second order. In 
the spin sector at weak coupling the backward scattering 
term is marginally irrelevant if C/ > 2V^ so that the spin 
excitations are gapless. For U < 2V this term becomes 
marginally relevant and a spin gap As appears. However, 
this gap only opens up exponentially slow, i.e., for a fixed 
U and V >Vc we expect 



As 



VV"- Kexp{-const/(F - VJ} 



(7) 



with Vc ~ U/2 at weak coupling. The phase transition 
in the spin sector is therefore of Kosterlitz-Thouless (KT) 
type. 

As Nakamura [6] first noted, there is no symmetry 
which fixes the amplitude of the Umklapp and backward 
scattering terms to be the same. So although these am- 
plitudes are identical to first order in the interaction pa- 
rameters, one would expect in general that they start 
to differ once higher order corrections are taken into ac- 
count. In this case an additional phase between the SDW 
and CDW phases would occur. As already outlined in 
the introduction different methods have given strong ev- 
idence that an additional phase with BOW order does 
indeed exist although some controversy about the extent 
of this additional phase remains. ^'^"^^ 

In the following, we will first develop a criterion which 
allows us to determine the second order line where the 
charge gap closes with high precision from thermody- 
namic data. Next, we will consider the KT-type tran- 
sition where the spin gap opens. Finally, we will provide 
some direct evidence that the new phase has long-range 
BOW order at zero temperature and does not extend 
beyond the tricritical point. 

To determine the line in the U, V-phase diagram where 
the charge gap closes, we consider the charge susceptibil- 
ity Xc- If a charge gap Ac exists, Xc temperatures 
T ^ Ac is described by Eq. (3). In the low temperature 
regime, Xc therefore will be larger the smaller the charge 
gap is. According to Eq. (6) we therefore expect the fol- 
lowing behavior of XciTo,V) at fixed U and fixed low 
temperature Tq: UV <Vc then Xc{To, V) increases with 
increasing V whereas Xc(Tq, V) decreases with increasing 
V ifV> Vc. 

For high temperatures Tq 3> 1, on the other hand, 
Xc{Tq,V) will always decrease with increasing V as can 
be easily seen from a high temperature expansion. Up to 
second order in 1/T we find 



XciT » 1) 



1 

2T 
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1 

2T 



(U/2 + V) 



(8) 



For V < Vc we therefore have the situation that 
dxc/dV > for T < 1 and dxc/dV < for T > 1. The 
compressibility curves for different V < Vc at fixed U 
therefore have to cross at least at one point. For V > Vc, 
on the other hand, we have dxc/ dV < for high as well 
as for low temperatures so that no crossing is expected. 
The different behavior of the compressibility curves for 
V > Vc and V < Vc is a very efficient criterion to deter- 
mine Vc as is shown in Figs. 11, 12 for the cases U = 2 
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FIG. 11: (Color online) Charge susceptibility Xc for U = 2. 
Upper panel: Xc for V = 0.5, 0.6, • ■ • ,1.0 (in arrow di- 
rection). The inset shows a zoom of the region around 
the crossing point at T* ~ 0.12. Lower panel: Xc for 

V = 1.04, 1.1, 1.2, • ■ ■ , 1.6 (main, in arrow direction) and 

V = 1.04, 1.06, 1.08, 1.1, 1.12, 1.13, 1.2, 1.3 (inset, in arrow di- 
rection) . 
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FIG. 12: (Color online) Charge susceptibility Xc for U ~ 
4. Upper panel: Xc for V = 1.5, 1.6, • • • ,2.1,2.15,2.17 (in 
arrow direction) . The inset shows a zoom of the region around 
the crossing point at T* ~ 0.54. Lower panel: Xc for V — 
2.16, 2.17, 2.2, 2.3, 2.4, 2.5 (main, in arrow direction) and V = 
2.16, 2.17, 2.18, 2.19, 2.2 (inset, in arrow direction). 



and [/ = 4, respectively. From Fig. 11, lower panel, we 
see that the first curve crossing the ones for larger V is 
the one for V = 1.04. We therefore find K = 1.05 ± 0.01 
for U = 2. In principle, the critical point can be deter- 
mined with this method even more accurately. Similarly 
for U = 4, the inset of the lower panel of Fig. 12 shows 
that the first curve crossing is F = 2.16 which leads to 
the estimate Vc = 2.165 ± 0.005 in this case. Both crit- 
ical values are in good agreement with the most recent 
zero temperature DMRG calculation.^^ Another interest- 
ing point is that in both cases the curves for V < Vc do 
not only cross but do so at one well defined point. I.e., 
there is a well defined temperature T* where dxc/dV « 
for all V. Similar well defined crossing points have also 
been observed in other systems and other thermodynamic 
quantities, as for example, the specific heat.""^^'^" 

For the spin susceptibility, Xs , there is only a spin gap 



FIG. 13: (Color online) Spin susceptibility ai U — 4 for differ- 
ent V = 1.8, 1.9,2.0,2.02, ••• ,2.1,2.15,2.2,2.3, ••• ,2.7. In- 
set: Spin gap As at ?7 = 4 as a function of V. The gap is well 
fitted by ~ 2.52^1^ - 2.02 • exp[-0.41/(y - 2.02)]. 



above the KT transition. For V > the temperature 
dependence of the spin susceptibility is then again given 
by Eq. (3). dxs/dV < for all temperatures so that 
the curves do not cross. The same is true for V < V^^: 
In the low-temperature limit the spin sector is then de- 
scribed by conformal field theory and 



(9) 



The spin velocity increases with increasing interaction 
strength so that again dxs/dV < for all temperatures. 
Therefore no qualitative change happens at the transition 
line. In principle, one can try to use the fact that there 
is universal scaling of certain ratios of thermodynamic 
quantities in the conformal regime. The entropy is given 
by Eq. (5) so that 



lim — — 

T^O TXs 



2^ 
3 



(10) 



is universal in the regime with gapless spin excitations. 
However, these formulas are only valid at temperatures 
T ^ Ac. Because the spin gap opens close to the point 
where the charge gap vanishes, this criterion turns out to 
be useless for our numerical calculations. We therefore 
have to determine the KT line by directly extracting the 
gap from the susceptibility curves. As an example, we 
consider again the case U = 4 (see Fig. 13). For small 
V the behavior is qualitatively consistent with Eq. (9) 
whereas for large V a spin gap is clearly visible. Fitting 
the low temperature part of the curves where a gap is 
present using Eq. (3) we can extract As as a function 
of V as shown in the inset of Fig. 13. Here, the error 
bars are obtained by varying the fit region. Another fit 
according to Eq. (7) then yields V/^ « 2.02 ±0.06 where 
the error estimate stems again from a variation of the fit 
region. Within the estimated errors we therefore obtain 
strong evidence that Vc ^ V^"^ for U = 4, i.e., that 
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FIG. 14; (Color online) Phase diagram as obtained by TMRG. 
The dashed line denotes U — 2V. The upper line describes 
the phase boundary of the CDW phase. The related error is 
always smaller than the symbol size. The error of the KT 
phase transition (lower line) is shown. 



that the transition into the CDW phase occurs. From 
field theoretical considerations it is then expected that 
the phase enclosed by these two transition lines is a Mott 
state with some dimerization, also called a bond-order 
wave (BOW) state. Because such a dimerization does not 
break any continuous symmetry, true long-range order at 
zero temperature will occur even in one dimension. This 
means that for the correlation function 



F{r) ^ i-iy {{AoAr) - {Arf) 



(11) 



with Ar — 5'^5'^+i or Ar = J2ai'^l,a'^r+i,cr + h.c.) wc havc 
linv^oo F{r) = const ^ 0. With' the TMRG algorithm 
there are different possibilities to detect this order. First, 
next-leading eigenvalues of the QTM allow it to calculate 
correlation lengths easily. In an asymptotic expansion 
of a two-point correlation function with operator Or we 
obtain 



{OiOr) - {Oi){Or) =J2M^e-^ 



(12) 



we have two separate phase transitions. Following this 
procedure to determine the second order and the KT 
transition lines for other values of [/, we obtain the phase 
diagram discussed in the next subsection. 



D. Phase diagram 

Our phase diagram, shown in Fig. 14, is very similar 
to the one obtained in the most recent zero-temperature 
DMRG study.^^ There is a first order transition line for 
{U,V) values above the tricritical point {Ut,Vt), sepa- 
rating the SDW and CDW phases. Below the tricriti- 
cal point we have a KT-type transition line where the 
spin gap opens and a second order phase transition line 
where the charge gap disappears. The nature of the so 
called BOW phase enclosed by the two transition lines 
is discussed in more detail in the next subsection. There 
is some quantitative difference between our study and 
Zhang's DMRG study in the location of the tricritical 
point though. We find Ut = 6.7 ± 0.2, Vt = 3.5 ± 0.1, 
whereas he found Ut « 7.2, Vt « 3.746. Both values are 
considerably larger than the ones found in QMC calcu- 
lations, Ut = 4.7 ± 0.1, Vt = 2.51 ± 0.04 in Ref. 7, and 
Ut G [5, 5.5] in Ref. 8. We also note that our phase dia- 
gram disagrees with that obtained in an earlier DMRG 
calculation^*^ where the BOW phase was restricted to the 
first order phase transition line (SDW-CDW) extend- 
ing from above the tricritical point estimated to be at 
Ut = 3.7 ±0.2 up to C/ w 8. 



E. Long-range BOW order and extent of the BOW 

phase 

From the phase diagram. Fig. 14, we see that the spin 
gap opens starting from the SDW phase and only after 



with correlation lengths ^q, and wave vectors ka given by 
Ao 



= In 



^.^arg('^ 



(13) 



where Aq is the largest eigenvalue of the QTM and Aq, an- 
other eigenvalue. A correlation length obtained accord- 
ing to Eq. (13) will show up in the asymptotic expan- 
sion (12) if the corresponding matrix-element Ma, which 
can also be calculated with the TMRG algorithm, ^^'^^ 
is nonzero. In the long distance limit, the behavior of 
the correlation function will be determined by the largest 
correlation length ^ with nonzero matrix-element. 

If the correlation function decays algebraically, this 
correlation length will diverge like ^ ~ 1/T. If, on the 
other hand, the correlation function decays exponentially 
even at zero temperature then ^ stays finite. Finally, for 
a correlation function showing true long-range order at 
zero temperature the correlation length will diverge like 



exp(A/r) 



T 



(14) 



where A is the gap for the corresponding excitations. 

In Fig. 15(a) we show, as an example, the leading 
SDW, CDW and BOW correlation lengths for [/ = 6 and 
V = 3.16. Here the leading SDW and CDW correlation 
lengths stay finite whereas the BOW correlation length 
diverges faster than 1/T indicating long-range BOW or- 
der at zero temperature. 

Another possibility to detect the BOW order with the 
TMRG algorithm is to calculate static susceptibilities 



Jo 



dT(Oo{0)Or{T)) (15) 



again for some operator Or ■ For true long-range order the 
corresponding xil) "^iH diverge exponentially with tem- 
perature, whereas xil) "^iH to a constant (zero if the 
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FIG. 15; (Color online) (a): Leading SDW, CDW and BOW 
correlation lengths plotted as for U — 6 and V — 3.16. 
(b)-(d): alternating static susceptibilities for the longitudinal 
spin, the density, and the kinetic energy k = X^o.('^r,<TCr+i,o- + 
h.c.) for U — 6 and different V, respectively. 



operator is conserved) for short-range order. The situa- 
tion is, however, complicated if the correlation function 
shows quasi long-range order, i.e., decays algebraically. 
Here we want to consider the case that only one sort of 
excitations is gapless, say the spin excitations. From con- 
formal field theory it is known that the corresponding 
algebraically decaying correlation function in the long- 
distance limit r 3> 1 will behave as 



(Oo(O)a(T)) ^ CXp 



-27rT.i; 



(16) 



X exp[— i/cr] exp[27rri((i+ — d )t] . 

Here v is the velocity of the elementary excitations, x ~ 
(i+ + the scaling dimension, the conformal weights, 
and k the characteristic wave vector. The r-integral for 
the static susceptibility xoo{k) can then be calculated 
explicitly and is given by 



/-/3 

/ dTexp[27rri(d+ - d-)T] 
Jo 



g27ri(d+-d ) _ 

2mT{d+ - d- ) 



(17) 



If the conformal spin s — c?+ — d~ is a non-zero integer - 
this is the case for any type of particle-hole excitation - 
the integral is zero and this part of the correlation func- 
tion does not contribute. If, on the other hand, s = 
then there is no time dependence in (16) and the inte- 
gral (17) yields just 1/T. The static susceptibility in 
the case of zero conformal spin will therefore scale as 
Xoo{k) ~ x'2a;-2 jj-^ particular, for the alternating part 
of the longitudinal spin-spin correlation function we have 



d+ = d~ = 1/4 leading to xs^S^M ~ ^/T. Note, how- 
ever, that for a; > 1 the long-distance asymptotics is no 
longer sufficient to discuss the behavior for T 0. In 
this case xoo{k) const for a non-conserved operator 
in general as in the case of exponentially decaying corre- 
lation functions discussed above. 

In Figs. 15(b), 15(c), 15(d) we show for J7 = 6 and 
different V the alternating static susceptibilities for the 
longitudinal spin, the density, and the kinetic energy, re- 
spectively. From Fig. 15(b) we conclude that a spin gap 
develops for V > 3.1. However, for V = 3.1 and V — 3.16 
there is still no long-range charge order (see Fig. 15(c)), 
i.e., an intermediate phase does exist. In Fig. 15(d) 
we see that at least for V = 3.16 this phase has long- 
range BOW order which is consistent with the correlation 
lengths shown in Fig. 15(a). Fitting the BOW correla- 
tion length using Eq. (14) we extract a rather small dimer 
gap A ~ 0.08. For fixed U the dimer gap is expected to 
decrease with decreasing V so that possible long-range 
bond order is detected most easily close to the transition 
into the CDW phase. In Fig. 16 the leading BOW corre- 
lation lengths for several U, values just below this tran- 
sition line are shown. For {U,V) = (4.0, 2.14), (5.5, 2.9), 
and (6.0,3.19) the correlation lengths diverge exponen- 
tially and we obtain the dimer gaps A « 0.01, 0.03, and 
0.08, respectively. As expected, A decreases with de- 
creasing U making it difficult to show the exponential 
divergence of the BOW correlation length for U < A be- 
cause temperatures below T ~ 10~^ are not easily ac- 
cessible by the TMRG method. Nevertheless, it is clear 
that the whole phase enclosed by the two transition lines 
shown in Fig. 14 must have long-range BOW order. For 
{U,V) = (7.5,3.9), on the other hand, we would expect 
A > 0.1 if BOW order does exist as found in Ref. 8 so 
that an exponential divergence should already become 
obvious at T ~ 0.1. However, down to T « 0.02 we see 
no indication of such a behavior, instead the BOW corre- 
lation length seems to diverge exactly as 1/T indicating 
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that wc arc in the SDW phase. This is supported by 
the data in the inset of Fig. 16 showing that the leading 
SDW correlation length also diverges like 1/T whereas 
the CDW correlation length stays finite for T ^ 0. In- 
terestingly, the BOW correlation length is larger than 
the SDW correlation length. We also confirmed that for 
[U,V) = (7.5,3.92) we are already in the CDW phase. 
We therefore conclude that for U = 7.5 no BOW phase 
exists. Instead, a direct first order transition from the 
SDW to the CDW phase occurs. 

IV. SUMMARY AND CONCLUSIONS 

Wc studied the thermodynamics of the half-filled one- 
dimensional extended Hubbard model using a TMRG al- 
gorithm. The focus was put on identifying the various 
phase transitions by considering thermodynamic quan- 
tities which are usually easy to measure like the uni- 
form magnetic susceptibility, Xs, or the isothermal com- 
pressibility, Xc- For strong coupling we calculated the 
charge gap in the SDW as well as charge and spin gap 
in the CDW phase in lowest order perturbation theory. 
The theoretical results were confirmed by TMRG calcu- 
lations of Xs and Xc- In the weak coupling regime where 
the phase transitions are continuous we showed that x,,- 
curves for a fixed U and different V as a function of 
temperature cross in one well defined point if measured 
in the SDW or BOW phase. In the CDW phase, on the 
other hand, no crossing occurs. We used this criterion 
to determine the boundary of the CDW phase with high 
accuracy. The KT transition line, on the other hand, 
where the spin gap starts to open exponentially slowly 
is very difficult to determine from thermodynamic data. 



Universal scaling relations obtained from conformal field 
theory for the magnetic susceptibility and the specific 
heat in the SDW phase turned out to be useless for this 
purpose. These scaling relations are only valid at tem- 
peratures T <c Ac which are not accessible by TMRG 
because the charge gap Ac is already very small near 
the KT transition. We could, however, show that ex- 
tracting the spin gap from the magnetic susceptibility 
where it is large enough and fitting it to a field theory 
formula does allow to determine the transition line rea- 
sonably well. In particular, the results clearly confirm 
that the two transition lines do not coincide and that 
an intermediate phase exists. By studying correlation 
lengths and static susceptibilities we confirmed that this 
additional phase has long-range bond order. We were 
also able to determine the tricritical point accurately and 
found Ut = 6.7 ± 0.2, Vt = 3.5 ± 0.1. Furthermore, we 
showed that above this point long-range bond order does 
not exist. Instead we find that BOW correlations can be 
dominant in this regime while still decaying algebraically 
at zero temperature. The resulting phase diagram is in 
good quantitative agreement with the most recent zero 
temperature DMRG study. '^-'^ However, it does not agree 
with the phase diagram found in Ref. 8 with the BOW 
phase existing even above the tricritical point. 
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